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(57) Abstract: An apparatus for estimating a plurality of channels from a plurality of distincttransmitting points, wherein an input 
signal includes * superposition of signals based on a first pilot sequence LransrniltaYe I'rtsrn a first transmitting point and ona second 
pilot sequence transmiltable from a second iiansmilling point, the firstpilol sequence and the second pilot sequence having a phase 
shift to each odier.cornpriscs a transformer (103; 303), which is operative to translate phase shifts between the first and the second 
pilot sequence into a delay, so that a first group of a set of discrete values of the transformed signal includes channel information 
of a first channel and that a second group of a. set of discrete values includes channel information of a second channel, a. selector 
(10.5: 305) for selecting the first group from tnc set of discrete values and lor selecting the second group from the set of discrete 
values to obtain selected groups. End means (107; 315) for determining channel information for die plurality of channels based on 
the selected groups. The seieerion of groups at the output of the transformer for obtaining channel information for the channel allows 
the construction of a simplified receiver. 
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Apparatus and Method for Estimating a Plurality of Channels 

Description 

The present invention is in the field of telecommunications 
and, in particular, in the field of channel estimation in a 
multiple input scenario , in which a receiver receives 

signals from more than one transmitting antennas. 

The steadily-increasing demand for high data rates 
necessary for today's and future mobile radio applications 
require high data rate techniques efficiently exploiting 
the available band width or, in other words , the achievable 
channel capacity. Therefore, multiple input multiple 
output (MIMO) transmission systems have achieved 
considerable importance in recent years. MIMO systems 
employ a plurality of transmitting points, each of the 
transmit ring points having a transmit antenna, and a 
plurality of receiving points , each of the receiving points 
having a. receiving antenna, to receive signals being 
transmitted by the multiple transmitting points through 
different communication channels . In MIMO techniques , 
where the signals impinging from several transmitter 
antennas need to be separated, space - time codes or 
special multiplexing methods are used. 

The signals impinging on each receive antenna are the 
s upe r - po s i t i on of the signals from N T antennas, where N T 
denotes a number of transmitting points . This implies new 
challenges for channel estimation. Channel parameters, like 
a channel impulse response or a channel transfer function 
are required for subsequent processing of the received 
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data. While the separation of the signals corresponding to 
several transmitting points, each of them having a transmit 
antenna, is a challenging task, the extension from a 
receiver having one antenna to a system with several 
receive antennas is straight forward, as long as the 
signals are mutually unccrrelated . The structure of the 
channel estimation units is independent of the number of 
receive antennas N R . The extension from a multiple input 
single output (MISO) system to a MI MO system is to employ 
N R parallel channel estimation units, one for each 
receiving point (receive antenna) . 

The use of coherent transmission technicues in wireless 
systems requires estimation and tracking of the mobile 
radio channel. Since the signals transmitted, from multiple 
transmit antennas are observed as mutual interference , 
channel estimation for MIMO systems is different from, the 
single transmit antenna scenario. MIMO systems can be used 
with a multi-carrier modulation scheme to further improve 
the communication capacity and quality of mobile radio 
systems. A prominent representative of multi-carrier 
modulation techniques is the orthogonal frequency division 
multiplexing {OFDM) technique . A block diagram of an OFDM 
system is shown in Fig. 7. 

Multi carrier modulation in particular orthogonal frequency 
division multiplexing (OFDM) has been successfully applied 
to a wide variety of digital communication systems over the 
past several years . In particular for the transmission of 
large data rates in a broadcasting scenario (e.g. digital 
TV) , OFDM' s superior performance in transmission over 
dispersive channels is a major advantage. OFDM has been 
chosen for various digital broadcasting standards, e.g. DAB 
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5 or DVB-T. Another wireless application of OFDM is in high 
speed wireless local area networks (WLAN) . 

OFDM was first introduced in the 1960s. An efficient 
demodulation utilising the discrete Fourier transform (DFT) 

10 was suggested by. S . Weinstein and P. Ebert, "Data 

Transmission by Frequency Division Multiplexing Using the 
Discrete Fourier Transform", IEEE Transactions on 
Communication Technology,, vol. COM-19, pp. 628-634, Oct. 
1971. By inserting a cyclic prefix into the guard interval 

15 (GI) longer than the maximum delay of the channel, inter- 
symbol interference (ISX) can be eliminated completely and 
the orthogonality of the received signal is preserved. 
Since future mobile communication systems should support 
data rates several times higher than current systems, 

20 multi-carrier systems with proper coding arid interleaving 

offer both efficient implementation through the application 
of the Fast Fourier Transform { FFT} and sufficient 
robustness to radio channel impairments. 

25 Another OFDM-based approach, termed multi-carrier code 
division multiplex access (MC-CDMA) , were spreading in 
frequency direction as has been introduced in addition to 
the OFDM modulation, as described in K. Fazel and L. Papke, 
"On the Performance of Con vo 1 u t i o na 1 ly - Coded CDMA/OFDM for 

30 Mobile Communication Systems", in Proc. IEEE Int. Symposium 
on Personal, Indoor and Mobile Radio Communications 
(PIMRC'93), Yokohama, Japan, pp. 468-472, Sep. 1993. MC- 
CDMA has been, deemed to be a promising candidate for the 
downlink of fourth generation systems . Moreover, a MC/CDMA 

35 system with a variable spreading factor has been proposed 
as described in H . Atarashi and M. Sawahashi, "Variable 
Spreading Factor Orthogonal Frequency and Code Division 
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Multiplexing (VSF-OFCDM) in 3 rd International Workshop on 
Multi-Carrier Spread-Spectrum & Related Topics (MOSS 
2001), Obe rp fa f f e nhc fen, Germany, ■ Sep. 2001. 

For OF DM -based MIMO systems, one OFDM modulator is employed 
on each transmitting point, while OFDM demodulation is 
performed independently for each receiving point. For 
OFDM, the signal stream is divided into N c parallel sub- 
streams. The i tri sub-stream commonly termed i tn sub-carrier 
of the 1 th symbol block {OFDM symbol) is denoted by X Iri . 
After serial to parallel conversion (S/P) performed by a 
S/P converter 701 an inverse discrete Fourier transform 

!DFT} with N FTT points is performed by an IFFT transformer 
7 03 on each block and, subsequently, the guard interval 

(SI) having N GI samples is inserted by a GI block 7 0S to 
obtain a signal x i>n after parallel to serial (P/S) 
conversion performed by a P/S converter 7 03. After digital 
to analogue (D/A) conversion, the signal x(t) is 
transmitted over a mobile radio channel with an impulse 
response h ( t , . t ) . The received signal at receive antenna v 
consists of superimposed signals from. transmitting 
points. Assuming perfect synchronisation, the received 
signal impinging at receive antenna v at sampling instants 
t=[n+lN syKi ] T 3pl is obtained. 

where n (t) represents additive white Gaussian noise, and 
Nsym=N FFT +K GI accounts for the number of samples per OFDM 
symbol . The * operator denotes convolution. The signal 
yi,n received by the receiver is first serial to parallel 
(S/P) converted by a S/P converter 709 and the guard 
interval is removed by a GI block 711. The information is 
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5 recovered by performing a discrete Fourier transform (DFT) 
on the received, block of signal samples (in Fig. ? a FFT 
transformer 713 is used) to obtain the output of the OFDM 
demodulation Yj., r in the frequency domain. The received 
signal at receive antenna v after OFDM demodulation is 
10 given by 



15 symbol and the channel transfer function (CTF) of transmit 
antenna u, at sub-carrier i of the 1 th OFDM symbol, 
respectively. The term N 1(i accounts for additive white 
Gaussian noise (A5?GN) with zero mean and. variance N c . 

2 0 When transmitting an OFDM signal over a multi-path fading 
channel, the received signal will have unknown amplitude 
and phase variations. For coherent transmission, these 
amplitude and phase variations need to be estimated by a 
channel estimator. 

.25 

In the following, reference is made to pilot symbol-aided 
channel estimation (PACE), where a sub-set of the 
transmitted data is reserved for transmitting known 
information, termed "pilot symbols". These pilots are used 
30 as side information for channel estimation. 

To formally describe the problem, the received pilot of 
OFDM symbol ID, at the {iD ; } th sub-carrier 




where Xjf 



and 



d H^ v) denotes the t: 



.ransmitted information 
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where X^ jDj , and fljj^ denotes the transmitted pilot 
symbol and the channel transfer function (GIF) of transmit 
antenna u, at sub-carrier i = T.D f of the / - 7d* OFDM symbol, 

10 respectively. It is assumed chat the CTF varies in the 1 
and in the i variable, i.e. in time and in frequency. The 
term N jD j D ^ accounts for additive white Gaussion noise. 
Furthermore, 1 represents the number of OFDM symbols per 
frame, N c is the number of sub-carriers per OFDM symbol, D- 

15 and D t denote the pilot spacing in frequency and time, and 
% is the number of transmit antennas. The goal is to 
estimate for all {l F i,u} within the frame Y 1(i is 

measured. Additionaly, the symbols Xjf at the location 
(l.i)^(jD,jD f ) are known at the receiver. The channel 

20 estimation now includes several tasks: 

1. The separation of N T superimposed signals, 

2. Interpolation in case that D t or D f are larger than 
one, and 

25 3. Averaging over the noise N 7dJD/ by means of exploiting 



In order to estimate given Y 7DjDf f there are N c 

equations with N C N T unknowns, when one OFDM symbol is 
30 considered. Thus, a straight- forward solution of this 
linear equation system does, in general, not exist. By 
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transforming Y- D jD to the time domain, the runber of 

unknowns car. be reduced, making it possible to solve the 
resulting equation system in the time domain. This 
approach has the advantage that DIT-based interpolation, 
which is a standard technique, can be combined with 
estimation and separation of N? superimposed signals in one 
step, resulting in a computationally efficient estimator. 

For time domain channel estimation for MIMO-OFDM systems, 
the received pilots of one OFDM symbol Y- !D . fD are pre- 

multiplied by the complex conjugate of the transmitted 

pilot sequence X~ D - fD , for 1 < i < N' p . Then the result 

isttransformed into the time domain via an W' P -point IDFT . 
Subsequently, the N T superimposed signals are separated by 
a matrix, inversion. The time domain channel estimate is 
obtained by filtering the output of the IDFT operation with 
a finite impulse response (FIR) filter. The DFT-based 
interpolation is performed simply by adding N c -Q zeros for 
the channel impulse response (CIR) estimates f i.e. to 
extend the length of the estimate of length Q to N c 
samples. This technique is called of zero padding . An N ' p - 
point DFT transforms the CTR estimate of the pilots to the 
frequency response estimate of the entire OFDM symbol. 

Estimators based on discrete Fourier transform ( DFT) have 
the advantage that a computationally efficient transform in 
the form of the Fourier transform does exist and that DFT 
based interpolation is simple . 

The performance of the estimation in general, is dependent 
on the choice of the pilot symbols . It is desirable to 
chose a set: of pilot sequences, which minimises the minimum 
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5 mean squared error (MMSE) criterium (which is a measure of 
the performance) and the computational complexity of the 
estimator. Estimators based on the least squares {LS) and 
the MMSE critariurr. for OFDM-MIMO systems have been 
systematically derived by Y. Gong and K. Letaief in: "Low 
10 Rank Channel Estimation for Space-Time Coded Wideband OFDM 
Systems," i::. Proc. IEZE Vehicular Technology Conference 
(VTC f 2001-Fall) , Atlantic City, USA, pp. 722-776, 2001. 

I. Barhumi et al describe in: " Optimal Training sequences 

15 for channel estimation in MIMO OFDM systems immobile 
wireless channels", International Zurich Seminar on 
Broadband Communications (IZS02) , February, 2002 a channel 
estimation and tracking scheme for MIMO OFDM systems based 
on pilot tones. In particular, the authors describe a 

20 channel estimation scheme based on pilot tones being 

orthogonal and phase-shifted to each other. Although the- 
pilot symbols described in 7. he above-cited prior art allow 
an accurate channel estimation, an enormous computational 
complexity at the receiver is required in order to perform 

25 matrix inversions required by the channel estimation 

algorithm. Due to this high computational complexity, the 
estimation scheme described in the above prior art document 
cannot be implemented at low cost, so that the disclosed 
algorithm, may not be suitable for mass-market mobile 

30 receivers . 

Y. Li, et {"Simplified Channel Estimation for OFDM Systems 
with Multiple Transmit Antennas," IEEE Transactions on 
Wireless Communications, vol. 1 , pp. 67-75, Jan. 2002) , 
35 proposed a channel estimation scheme for OFDM with multiple 
transmit antennas which is based on the DFT transform. In 
particular, Li discloses a method for generating pilot 
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symbols to be transmitted by multiple transmit and receive 
antennas and to be exploited at the receiver for channel 
estimation. These pilot symbols are generated by 
multiplying a training sequence having good timing and 
frequency synchronisation properties by a complex sequence 
introducing an additional phase shift between the pilot 
symbols and between the subsequent values of each pilot 
symbol, as well. To be more specific, each value of a 
training sequence is multiplied by a complex factor, which 
introduces a phase shift , wherein the phase shift is 
dependent of a number being assigned to the value being 
multiplied, on a number assigned to the corresponding- 
transmitting point and. a total number of transmitting 
points. The pilot symbols are orthogonal and phase shifted 
to each other. The pilot symbols are modulated by an OFDM 
scheme and transmitted through a plurality of communication 
channels. At a receiver,, which is one cf a plurality of 
receivers, a signal being received includes a super- 
position of the plurality of transmitted signals through 
the plurality of communication channels, Li et al presented 
further a design rule for the pilot tones based on phase- 
shifted sequences which is optimum in the mean squared 
error (M3E) sense. Moreover, a matrix inversion, which .i s, 
in general, required for the estimator, can be avoided by 
choosing orthogonal pilot sequences. However, due to a 
difficulty of obtaining perfect orthogonality between 
training sequences, matrix inversions may be necessary. 

Fig . 9 shows prior art channel estimation scheme as t ought 
by Li, where the case of two transmitting antennas is 
considered. 
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The prior art channel estimator includes a plurality of 
multipliers, wherein Fig. 9 shows only three multipliers 
being associated with the V. th value of a n th received 
sequence r [n, k] . A first multiplier 301, a second 
multiplier 903 and a third multiplier 905 arranged in 
parallel include, first and second inputs and outputs, 
respectively. The output of the first multiplier 901 is 
connected to a first inverse fast Fourier transform ( I FFT } 
block 907, the output of the second multiplier 903 is 
connected to a second IFF? block 909 and the output of the 
third multiplier 305 is connected to the third I FFT block 
911 . It should be mentioned here that in total, K 
multipliers are connected to each I FFT block, wherein K 
denotes a length, of a received sequence in the frequency 
domain, and a total number of 3K input signals are provided 
to the three I FFT blocks. Each of the I FFT blocks 907, 909 
and 911 is operative to perform an inverse fast Fourier 
algorithm applied to K input values, respectively. 
Furthermore, each of the I FFT blocks 907 f 909 and 911 
includes a number of outputs, wherein only the first K 0 
outputs of each I FFT block are used.. The respective 
remaining outputs are, for example, connected to ground. 

K Q outputs of the first I FFT block 907 are connected to a 
first estimation block 913 and the first K 0 outputs of the 
third I FFT block 911 are connected to a second estimation 
block 915. The K c outputs of the second I FFT block 909 are 
connected to the first estimation block 913 and to the 
second estimation block 915, respectively. The first 
estimation block 913 and the second estimation block 915 
have K 0 outputs, each of the outputs being connected to a 
corresponding filter 917 of a plurality of filters, each of 
the filters having an output, respectively. The K 0 outputs 
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5 of the filters 317 corresponding to the first estimation 
block 913 are connected to a first Fourier transform (FFT) 
block 917 and the K 0 outputs, of. the filter 917 
corresponding to the second estimation block 915 are 
connected to a second FFT block 921 . The first FFT block 

10 919 and the second FFT block 921 have K outputs, where K 
is, as stated above, the number of sub-carriers. 
Furthermore , due to the simplified algorithm described by 
Li, the outputs of the first filters 917 corresponding to 
the first estimation block 913 are connected to the second 

15 estimation block 915 and the outputs of the filter 917 
corresponding to the second estimation block 915 are 
further connected to the first estimation block 913, so 
that a plurality of feedback loops is established. 

20 As stated above, Fig. 9 shows an example of the prior art 
estimator for the case of two transmit antennas, so that 
the received signal r[n,k] is a superposition of two 
transmitted signals being possibly corrupted by channel 
noise . The received signal is split into two received 

25 signals by a splitting means not shown in Fig. 9, The 
copies of the received signals are then multiplied by 
complex conjugated signals corresponding to the respective 
transmit antennas. Moreover, the pilot symbol transmitted 
by the first transmit antenna is multiplied by a complex 

30 conjugated version of the pilot symbol transmitted by the 
second antenna. More precisely, the K values of the 'first 
copy of the received signal are multiplied by K values of 
the complex conjugated version of the pilot symbol 
transmitted by the first antenna. The K values of the 

35 second version of the received signal is multiplied by K 
values of the complex conjugated version of the pilot 
symbol transmitted by the second transmit antenna. 
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Furthermore, the K values of the pilot symbol transmitted 
by the first antenna is multiplied by K complex conjugated 
values of the pilot symbol transmitted by the second 
transmit antenna in order to obtain intermediate values 
required by the subsequent channel estimation algorithm. 

As stated above, all multiplications are performed in 
parallel, so that the K results from the multipliers 901 
are fed to the first IFFT block 907. K results from the K 
multipliers 903 are fed to the second IFFT block 909. K 
results from the K multipliers 905 are fed to the third 
IFFT block 911. Each respective IFFT block is operative to 
perform an inverse fast Fourier transform in order to 
transform the frequency domain input signals into time 
domain output signals. 

The first and the second estimation block, 913 and 915, are 
operative to perform a channel estimation algorithm based 
on the plurality of the input signals. More precisely, the 
first estimation block S13 receives 3Kc input signals to 
generate K 0 output signals corresponding to the channel 
impulse response of the first channel from the first 
transmit antenna to the considered receive antenna . The 
second estimation block 915 receives, in an analogue way, 
3K(i input signals to generate K 0 output values 
corresponding tc the second communication channel from the 
second transmit antenna to the receive antenna. The 
respective K 0 output values are then filtered by filters 
917. 

As stated above , the respective output signals from the 
filters are fed back to the first and second channel 
estimation blocks 913 and 915, since the channel estimation 
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5 blocks 913 and 915 are operative to estimate the channel 
impulse response of the respective communication channels 
based on previously-calculated values and on current values 
obtained from the IFFT blocks. Each estimation block 
applies an estimation algorithm where matrix by vector 

10 multiplications .instead of matrix inversions are performed 
in order to calculate desired channel impulse responses . 
After filtering and zero padding to a length required by 
the following fast Fourier transform, a channel transfer 
function of the first and of the second communication 

15 channels are obtained. 

As stated above, Li avoids matrix inversions by introducing 
an iterative scheme where matrix by vector multiplications 
appear and by exploiting the orthogonality of the pilot 

20 symbols. However, in order to calculate two channel 

impulse responses corresponding to the two communication 
channels, three inverse fast Fourier transforms, a 
splitting means and 3K multipliers are required. Moreover, 
the channel estimation algorithm applied by Li has still a 

25 high complexity due to the required matrix by vector 
multiplications. Hence, with an increasing number of 
transmit antennas r the complexity of the complicated 
estimation scheme proposed by Li rapidly increases due to 
the high number of complex valued multiplications. In 

30 addition, the multiplication of the two pilot symbols 

followed by an inverse Fourier transform is necessary in 
order to provide a plurality of intermediate values 
required for channel estimation. Hence, the estimation 
blocks 913 and 915 cannot operate independently, so that 

35 additional timing and control operations are necessary. 
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It is an object of the present invention to provide an 
improved concept for channel, estimation with reduced 
complexity. 

This object is achieved by an apparatus for estimating a 
plurality of channels in accordance with claim. 1 or an 
apparatus for providing a plurality of different pilot 
sequences in accordance with claim 2 4 or an apparatus for 
providing a number of different pilot sequences in 
accordance with claim 25 or a method for estimating a 
plurality of channels in accordance with claim 2 6 or a 
method for providing a plurality of different pilot 
sequences in accordance w.ir.h claim 27 or a method for 
providing a number of different pilot sequences in 
accordance with claim 28 or a method for providing a number 
of cifferent pilot sequences in accordance with claim 2 9 
computer programme in accordance with claim. 30. 

The present invention is based on the finding that a 
channel estimation based on orthogonal and phase-shifted 
pilot symbols with respect to each other can efficiently be 
performed by translating the phase shifts into delays, 
wherein an input signal including a superposition of a 
plurality of signals based on the pilot sequences 
transmitted through a plurality of communication channels 
is transformed into a transformed signal, so that the phase 
shifts between the plurality of pilot symbols are 
translated into delays . Hence, a set of discrete values of 
the transformed signals can be divided into a plurality of 
groups, each group having a delay with respect to the 
preceding groups, wherein each group represents channel 
information of a channel. In particular, it has been found 
that a single transformer can be applied for transforming 
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the phase shifts into delays r since the delayed groups 
actually include channel information, i.e. a channel 
impulse response of a communication channel, when a phase 
shift to delay transform is performed. Therefore, by 
exploiting the distinct delays in order to select the 
groups, the channel information can be obtained directly 
from, the selected groups using only a single transformer. 

For the sake of simplicity, in the following, only one 
receiving point being representative for the plurality of 
receiving points, will be discussed. However, it should be 
pointed out here that the following discussion and results 
can be applied to a plurality of receiving points in a 
straight-forward way. Furthermore, in the following an OFDM 
transmission system will be considered. However, the 
present invention can also be applied to channel estimation 
in single carrier modulation, systems or in multichannel 
modulation techniques different from OFDM. To this end, 
preferably a guard interval is used so that the transmitted 
single carrier symbol is prevented from being affected by 
received signal from inter-symbol interference. 

Referring again to the case of two transmitting points, an 
inventive apparatus for estimating a plurality of channels 
from a plurality of distinct transmitting points to a 
receiving point, derives an input signal from a signal 
received at the antenna of the receiving point . The input 
signal includes a superposition of a first signal based on 
a first pilot sequence transmit table from a first 
-ransmitting pcint and a second signal based on a second 
pilot sequence transmit table from a second transmitting 
point. The first pilot sequence and the second pilot 
sequence are orthogonal to each other within a 
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predetermined orthogonality range. If the first pilot 
sequence and the second pilot sequence are perfectly 
orthogonal to each other, then, an- inner product of these 
two sequences equals 0. However, since a perfect 
orthogonality is difficult to achieve, it is allowable that 
an absolute value of an inner product of the first pilot, 
sequence and of the second pilot sequence is greater than 
or equal, to 0 but preferably smaller than 0.2. Hence, the 
predetermined orthogonality range is defined by an interval 
having a first value equal to 0 and the last value equal to 
0.2. 

If the first and the second pilot sequences have an 
additional phase shift to each other - then a transformer 
for transforming the input signal into a transformed signal 
translates the phase shifts into a delay r so that a first 
group of a set of discrete values of the transformed signal 
includes channel information of a first channel from the 
first transmitting point to the receiving point and a. 
second group of the set of discrete values includes channel 
information of a second channel from the second 
transmitting point to the receiving point. Hence, as 
stated above, the first channel information is extracted 
simply by selecting the first group from the set of 
discrete values and the second channel information is 
extracted by selecting the second group from the set of 
discrete values . In order to obtain selected groups, the 
inventive apparatus further comprises a selector for 
selecting the first and the second group based on the delay 
between the first and the second group. Based on the 
selected groups, the channel information can be determined 
by a means for determining channel information. 
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5 In accordance with the present invent lor., an input signal 
is transformed by the transformer to translate, the phase 
shift into a delay. Preferably, the transformer is 
operative to perform a Fourier transform, since a phase 
shift in the frequency domain is translated to a delay in 

10 the time domain .by, for example, an inverse discrete 
Fourier transform and vice-versa. So if the pilot 
sequences are phase-shifted in the frequency domain, the 
resulting time domain signal after the inverse Fourier 
transform will be delayed. Furthermore, by assigning a 

15 certain phase shift; to each pilot sequence, the 

corresponding delay in the time domain will be dependent on 
the transmit antenna. Provided -hat N' p >N T Q f wherein 
denotes a number of discrete values of the input signal, 
and by using appropriate seta of phase shifted sequences, 

20 uhe N T superimposed signals can be simply separated by 
performing an inverse discrete Fourier transform. 

While the estimator front and disclosed by Li et al 
requires pre-multiplications and a plurality of ID FT 

2 5 operations, only one ID FT operation is required to separate 
Mt superimposed signals by the inventive estimation scheme. 
The remaining part of the receiver will not be affected by 
this modification. Although certain phase-shifted 
sequences are perfectly suited in time-domain for OFDM 

30 channel estimation, other orthogonal sequences, such a 

sequences derived from Hadamard matrices, may be also be 
applied to the proposed simplified receiver structure, with 
some modifications resulting in a modest additional 
processing in time domain, which will be described later. 



In general, a significant improvement of the proposed 
receiver structure compared to the prior art receiver 
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structure is the simplification, of the estimator front-end. 
This simplification has no negative effects on the 
performance of the channel estimator . Moreover, if the 
phase shifted pilot sequences are appropriately designed, 
then the channel estimation scheme can further be 
simplified, since phase shifted pilot sequences are a 
perfect match to a Fourier transform in terms of separating 
the N T superimposed signals without any substantial further 
processing . 

In addition, an increasing number of transmitting points 
has only a slight impact on a complexity of the inventive 
channel estimation scheme, since irrespective of the number 
of transmitting points, only one transformer is required to 
trans form the input signal into transformed signal so that 
the phase shifts are translated into delays. Accordingly,, 
only one selector is applied to select the selected groups. 
This is a significant complexity reduction compared to the 
prior art discussed above, since the complexity of the 
prior art channel estimation scheme rapidly increases with 
an increasing number of transmitting points due tc- the 
required, plurality of IDFT transformers in order to 
transform the plurality of (premultiplied) copies of the 
received signal, and in order to obtain a. further plurality 
of intermediate values required by the plurality of channel 
estimation blocks. 

Furthermore, the inventive apparatus does not comprise any 
splitting means since only one input signal is directly 
transformed which leads to a further complexity reduction. 
Moreover, since only one transformer is used, an optional 
preprocessing of the input signal prior to transformation 
(windowing, zero padding etc.) is simpler. 
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Furthermore, due to the inventive channel estimation 
scheme, significantly less multiplications or no 
multiplications at all have to be performed which leads to 
a further significant complexity reduction when the input 
signal is not pre-multiplied by a plurality of complex 
conjugated pilot sequences. If the pilot sequences are pre- 
multiplied by a scrambling sequence in order to achieve 
better transmission characteristics of the signal to be 
transmitted, then only the values of one input signal have 
to be pre-rr.ultipliec for removing the influence of the 
scrambling sequence. If nc pre multiplying is performed, 
then no multiplications prior to transforming are 
necessary. On the contrary, the prior art channel 
estimation scheme always requires a plurality of 
multiplications as described above. 

Furthermore, the inventive channel estimation scheme is 
simplified since the transformed signal contains actually 
the channel information. Therefore, no matrix inversions or 
matrix by vector multiplications are necessary which 
further reduces the complexity. 

Moreover, the inventive channel estimation scheme can be 
applied to any orthogonal sequences having different phase 
shifts to each other provided that these phase shifts are 
known at the transmitting point. Since these different 
phase shifts are transformed into delays, the inventive 
selector can always select proper groups from a plurality 
of delayed discrete values even when the discrete values 
belonging to a particular group are not succeeding each 
other. 
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5 Moreover, the performance of the proposed channel 

estimation scheme can be further increased by providing a 
plurality of different pilot sequences in accordance with a 
further embodiment of the present invention. In order to 
calculate an i th value of a pilot sequence to be 

10 transmitted by a. transmitting point from a number of 

transmitting points, the inventive apparatus for providing 
a plurality of different pilot sequences comprises a means 
for calculating a pilot sequence being operative to 
generate predetermined phase shifts between succeeding 

15 values of the pilot sequence to be transmitted and 

predetermined phase shifts between the pilot sequences as 
well. 

To be more specific, the phase shifts between succeeding 
20 values of the pilot sequence to be transmitted are 

generated, such that the phase shift between successive 
values of the pilot sequence is negative in a mathematical 
sense. Furthermore,, the phase shift between the succeeding 
values depends on a number associated with the 
25 corresponding value in the pilot sequence. 

In addition, in order to generate phase shifts between the 
plurality of pilot sequences, the phase shifts between 
succeeding values in a particular pilot sequence also 

30 depends on a number assigned to a transmitting point, the 
number being divided by a total number of transmitting 
points. Moreover, the generated pilot sequences are 
preferably scrambled by the same scrambling sequence, the 
length of the scrambling sequence being equal to the length 

35 of the pilot sequences, wherein each value of the pilot 

sequence to be generated is multiplied by a corresponding 
value of the scrambling sequence. Hence, all pilot 
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sequences include the same possibly complex valued 
scrambling sequence introducing an additional phase shift 
and amplitude modulation, Therefore, at the receiving 
point, the additional phase shift introduced by the 
scrambling sequence can be removed by applying a single 
multiplier by simply multiplying the input signal by a 
complex conjugated version of the scrambling sequence. 
Moreover, the inventive phase shift between each succeeding 
value of each pilot sequence is suited to an inverse 
discrete Fourier transform in order to translate the phase 
shifts into delays, since the inventive phase shift between 
succeeding values of each, pilot, sequence to be transmitted 
is negative in a mathematical sense and is hence a perfect 
counterpart of the IDFT (translating phase shifts into 
delays) since the XDFT introduces a phase shift in a 
positive mathematical sense. Hence, a simplified DFT based 
channel estimation scheme can be performed. 

Moreover, in order to simplify the channel estimation 
scheme used at the receiving point, suitable orthogonal 
sequences can be used as pilot sequences (pilot symbols) . 
In accordance with a further embodiment of the present 
invention, an apparatus for providing a number cf different 
pilot sequences to be used for channel estimation comprises 
means for generating a Hadamard matrix having a number of 
columns, each column being orthogonal to the other columns, 
wherein the number of columns corresponds to the number of 
different pilot symbols . Based on the columns of the 
Hadamard marrix, the different pilot sequences can be 
calculated by, for example, a cyclic extension of each 
column cf the Hadamard matrix, wherein each pilot, sequence 
is derived from another column of the Hadamard matrix. 
Since the generated pilot sequences have a predetermined 
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phase shift to each other, these phase shifts translated, 
into delays at the receiver by, for example, applying an 
inverse discrete Fourier transform to obtain discrete 
values and to select groups containing channel information 
from the delayed discrete values . Moreover, this approach 
is suited for multi-carrier cods division multiplex access 
transmission systems (MC-CDMA) , where spreading in 
frequency direction is introduced in addition to the OFDM 
modulation. In order to introduce spreading in frequency 
direction, orthogonal signals based on columns of a 
Hadaiuard matrix can be applied, which simplifies a channel 
estimation procedure, since the same sequence can be 
applied for spreading in the frequency direction and for 
reduced complexity channel estimation. The inventive 
approach is also suited for channel estimation in multiple 
access transmission systems, like FDMA (frequency division 
multiple access) or TDMA (time division multiple access) . 

Further embodiments of the present invention are described 
in detail with respect to the following figures, in which: 

Fig. 1 shows a block diagram of an inventive apparatus 
for estimating a plurality of channels; 

Fig. 2 demonstrates a preferred transformation of an 

input signal into a plurality of delayed groups; 

Fig. 3 shows a block diagram of an further inventive 

apparatus for providing a plurality of different 
pilot sequences ; 
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5 Fig. 4 shows a block diagram of a inventive apparatus 
for providing a number of different pilot 
sequences; 

Fig. 5 shows a block diagram of a further inventive 
10 apparatus for providing a number of different 

pilot sequences; 

Fig. 6 shows a pilot grid structure ; 

15 Fig. 7 demonstrates an OFDM modulation and demodulation 
scheme; and 

Fig. 8 shows a block diagram of prior art channel 
estimation schemes. 

20 

The apparatus shown in Fig. 1 comprises a multiplier 101 
having two inputs and one output . The output of the 
multiplier 101 is connected to a transformer 103 having a 
plurality of outputs for providing a transformed signa. The 

25 plurality of outputs of the transformer 103 is connected to 
a selector 105 living a control input 104 providing a delay 
information for proper group selection. The selector 105 
has a plurality of outputs, wherein a number of outputs is 
equal to a number of inputs of the selector 105 and, thus, 

30 equal to a number of outputs of the transformer 103. The 
plurality of outputs of the selector 105 is connected to a 
means 107 for determining channel information, the means 
107 having a plurality of outputs. 

35 The apparatus shown in Fig. 1 demonstrates the inventive 
low complexity channel estimation scheme based on pilot 

signals being orthogonal to each other within the 
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5 predetermined orthogonality range and having a phase shift: 
to each other. If the pilot sequences comprised of pilot 
symbols, each of the pilot sequences corresponding to one 
of a plurality of transmitting points, are scrambled by the 
same complex amplitude sequence, then the pilot sequences 

10 include a complex amplitude sequence and phase information. 
The transmitting points then transmit signals based on the 
pilot sequences through a plurality of communication 
channels, wherein prior to transmitting, the pilot: 
sequences are modulated using, for example, an OFDM 

15 modulation scheme. A signal received at a receiving point 
(not shown in Fig. 1) includes a superposition of the 
plurality of transmitted signals and possibly a plurality 
of superimposed distortion signals representing channel 
noise. After a demodulation performed by applying an OFDM 

20 demodulation scheme, a received signal shown in Fig. 1 is 
obtained. The multiplier 101 is, hence,- operative to 
multiply the received signal by a complex conjugated 
version of the amplitude sequence (scrambling sequence) in 
order to remove the amplitude information and phase shift 

25 introduced to each of the pilot sequences by scrambling. 
The received signal is a set of successive discrete 
received signal values and the multiplier 101 is operative 
to perform a single complex multiplication for each value 
of the set of successive discrete received signal values 

30 and a corresponding value of a set of successive discrete 

complex conjugated amplitude signal values, so that after a 
number of single complex multiplications, a set of 
successive discrete input values of an input signal is 
obtained, 

35 

The input signal includes a superposition of signals based 
on a plurality of pilot sequences, for example, based on a 
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first, a second and a third pilot sequence being- 
transmitted from a first, a second and a third transmitting 
point . Hence , the input signal contains phase shifts of 
the pilot sequences, wherein the phase shifts are not 
removed by further pre-multiplications . Instead, the input 
signal is provided to the transformer 103 being operative 
to translate the phase shifts into delays. 

The transformer 103 can be any transformer being able to 
translate phase shifts into delays. For example, the 
transformer 103 performs a Fourier transform or a fast 
Fourier transform or a discrete Fourier transform or an 
inverse Fourier transform or an inverse fast Fourier 
transform or an inverse discrete Fourier transform of the 
input signal. In particular,, the input signal is a set of 
successive input signal values and the transformer is 
operative tc transform the set of successive input signal 
values into the transformed signal, which is a set of 
discrete values, wherein a total number of discrete values 
of the set of discrete values of the transformed signal is 
equal to a total number of input signal values in the set 
of input signal values. 

Hence, the transformer is operative to perform a single 
transform algorithm applied to the set of discrete values, 
wherein the transform algorithm can be, for example, one of 
the algorithms mentioned above applied to the total number 
of input signal values. The length of the input signal 
corresponds to a length of a training or pilot sequence 
multiplied by a total number of transmitting points. 
Hence, the length of the transformed signal is equal to the 
length of the input signal. 
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The selector 105 selects a number of groups, the number of 
groups being equal to the total number of transmitting 
points in order to obtain selected groups, wherein each 
group consists of an equal number of successive discrete 
values of the set of discrete values of the transformed 
signal. The selector is operative to select, as a group, a 
number of successive discrete values, wherein the number of 
successive discrete values is equal to the total number of 
discrete values of the transformed signal divided by the 
total number of transmitting points, wherein the total 
number of discrete values is obtained by the transformer. 

For example, the selector is operative to select groups 
such that a first discrete value of the set of discrete 
values is a first discrete value in a first group at the 
output of the selector and a first discrete value of a 
second group of discrete values is a discrete value in the 
set of discrete values succeeding a discrete value having 
the highest order of the values selected to obtain the 
first group, wherein the first group and the second group 
have ccr re spending distinct delays, the delays being 
obtained by the transformer. If the input signal includes, 
in a superimposed manner, a third signal based on a third 
pilot sequence transmittable from a third transmitting 
point, the transformer translates a further phase shift of 
the third group into a further delay, so that a third group 
of the set of discrete values includes channel information 
of a third channel, from the third transmitting point to the 
receiving point and the selector selects the third group 
from the set of discrete values to obtain a selected third 
group having a further delay associated with the third 
group. Hence, if the input signal includes a plurality of 
superimposed signals corresponding to a plurality of 
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5 transmitting points, the selector 105 selects a further 
group of the plurality of groups on the basis of a delay 
corresponding to the further group as depicted in Fig. 1. 

The groups being selected by the selector based on. a delay 

10 information via the control input 104 or based on a fixed 
setting (no control input is required in this case) are 
provided to the means 107 for determining channel 
information based on the selected groups. If the pilot 
symbols are orthogonal to each other within the 

15 orthogonality range and if the pilot symbols have a 

predetermined phase shift with respect to each other, then 
each of the selected groups contains channel information 
corresponding to a particular channel from a particular 
transmitting point to the receiving point. In this case, 

20 the means 107 is operative to receive the selected groups 

and to output the selected groups, since each of the groups 
includes channel information of the particular channel, as 
will be shown later. In order to output the selected 
groups, the means 107 can be designed in such a way that 

25 the Inputs of the means 107 are directly connected to the 
outputs of the means 1C7. If there is no channel noise or 
If the channel noise Is negligible, then the means 107 for 
determining channel information outputs the selected 
groups, each group including, for example, a channel 

30 impulse response of a particular channel as channel 

information. Furthermore, the means 107 can include "a time 
frequency conversion means, for example, a DFT based 
transformer for obtaining a channel transfer function by 
transforming the determined channel impulse response of a 

35 particular channel into the frequency domain in order to 
obtain the channel transfer function as the channel 
information. 
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5 

Depending on pilot sequences applied to channel estimation, 
a group of the selected groups _ and a further group of the 
selected groups may include channel information of a 
channel and channel information of a further channel in a 

10 superimposed manner. Furthermore, depending on the design 
of pilot sequences, a group of rhe selected groups may 
include a superposition of the channel information of a 
channel being multiplied by a phase shift factor and of a 
further channel information of a further channel being 

15 multiplied by a further phase shift factor, wherein a 
further group of the selected groups may include a 
superposition of the channel information of the channel 
being multiplied by the complex conjugated phase shift 
factor and of the further channel information of the 

20 further channel being multiplied by the complex conjugated 
further phase shift factor, wherein the phase shift factor 
and the further shift factor are caused by the design of 
the pilot sequences. In this case, the means 107 is 
operative to determine the channel information of the 

25 channel and the further channel information of the further 
channel using the group and the further group by, for 
example, solving a set of linear equations and removing the 
additional phase shift factors. For example, the means 107 
applies phase shift operations and additions and/or 

30 subtractions to the group and to the further group in order 
to determine channel information. If, for example, pilot 
sequences based on columns of a Hadamard matrix are applied 
to channel estimation, then the phase shift factor is a 
function of a phase shift of 180° divided by the total 

35 number of transmitting points. Using this information, the 
known phase shifts are removed to determine the channel 
information. 
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5 

If the channel noise cannot be neglected, then the channel 
information in the selected groups is affected by a 
distortion signal related to channel noise. In order to 
suppress the distortion signal, the means 107 may comprise 

10 an estimator being operative to estimate the channel 

information from the corrupted channel information by, for 
example, applying a minimum mean squared error ( MMS E ) 
estimator , a least squares (LS) estimator, a maximum- 
likelihood (ML) estimator or a maximum a posteriori (MAP) 

15 estimator (being a particular embodiement of a ML 

estimator) . If the MMSE estimator is applied to estimate 
channel information, then the estimator performs filtering 
of the corrupted channel information, wherein the filter 
has filter coefficients, wherein each filter coefficient is 

20 determinable from, a multiplication of a first matrix by a 
second matrix, wherein the first matrix is a cross- 
correlation (In general; cross-covariance) matrix between 
channel information and corrupted channel information and 
the second, marrix is an inverse of an autocorrelation (in 

25 general; auto-covariance) matrix of the corrupted channel 
information. 

The pilot sequences applied to channel estimation in an 
OFDM transmission system are commonly used to modulate 

30 discrete carriers in a multi-carrier modulation scheme at a 
predetermined time instant corresponding to the time 
instant at -which pilot sequences are to be transmitted. 
Hence, the inventive apparatus may further comprise a 
controller for detecting the input signal in a received 

35 signal and for enabling rhe transformer, the selector or 
the means for detecting only when the input signal 
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5 including pilot sequences is detected in the received 
signal . 

Alternatively, the pilot sequences can be used to modulate 
the same pilot carrier in a multi-carrier modulation scheme 

10 in order to estimate time varying coefficients of a channel 
transfer function of a communication channel, wherein the 
coefficients of the channel transfer function are 
determined by the pilot carrier. Hence, the inventive 
apparatus can comprise a controller for detecting the pilot 

15 carrier in a received signal, including successive multi- 
carrier modulated symbols and for collecting values of the 
pilot carrier at time instants corresponding to the 
transmitting instants tc obrain the input signal, wherein 
the channel information obtained by processing the input 

20 signal using the transformer selector and the means for 

detecting includes a behaviour of a multi-carrier channel 
from transmitting points to receiving points defined by the 
pilot carrier over time. The channel information including 
the behaviour of the multi-carrier channels is f for 

25 example, related to a velocity of a mobile receiver 
causing, for example, a Doppler shift related to the 
relative velocity between transmitter and receiver. In 
order to estimate the Doppler shift, changes of a 
coefficient of the channel transfer function over time can 

30 be exploited. 

If the pilot sequences are transmitted, at. distinct time 
instances and are not succeeding each other, preferably an 
interpolation is performed in order to obtain, intermediate 
35 values of channel information between time instances 

corresponding to the time instantsat which pilot sequences 
are transmitted. Hence, the means for determining channel 
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information may comprise a further means for interpolation 
being operative to interpolate between a first channel 
information obtained at a first time instantsand a second 
channel information obtained at a second time instance, 
wherein the intermediate interpolated channel information 
can, for example, be determined by applying a well known 
Wiener interpolation, technique (Wiener filter) based on the 
first and the second channel information or, for example, 
by applying a Langrange interpolation scheme . 

Fig. 2 demonstrates a principle of exploiting the DFT 
properties by using phase-shifted sequences in the 
frequency domain , wherein the DFT algotithmus separates N T 
signals in the time domain, wherein W T corresponds to a 
number of transmitting points. The frequency domain iirout 
signal can include pilots having a spacing in frequency 
denoted by D ff as depicted in Fig. 2- The pilots are 
orthogonal to each other within the orthogonality range and 
have different phase shifts which can be translated into 
delays by applying the DFT-based transform yielding a time 
domain transformed signal. The transformed signal in the 
time domain includes N T groups including channel 
information like, for example, a channel impulse response 
of a particular communication channel . Each of the groups 
includes K P /N T discrete values, wherein N p denotes a total 
length of the input signal. The groups transformed by the 
transformer are time-delayed, wherein a second group has a 
delay with respect to a first group, so that the discrete 
values of the second group do not overlap with the discrete 
values of the first group. Accordingly, the N T th group 
corresponding to the N T transmitting point has the same 
delay with respect to a preceding group . N FFT denotes a 



WO 2IMI4,'HX<i7(>6 



32 



PCT/EP2003/00J211 



number of discrete FFT points being dependent on the pilot 
spacing in frequency domain. 

It is well known that a very efficient implementation of 
the DFT is the fast Fourier transform (FFT) algorithm. For 
optimum efficiency, the number of points of the FFT should 
be a power of 2, It is important to note that the 
interpolation ratio should be D f . 

Fig. 3 shows a conceptual block diagram of an inventive 
apparatus for estimating a plurality of channels from. a. 
plurality of distinct transmitting points to a receiving 
point . 

The apparatus shown in Fig. 3 comprises a frequency domain 
(FD) window 301 having N? inputs corresponding to N P 
discrete values of an input signal provided to the 
frequency window 301 and N P outputs connected to a l%~ I FFT 
transformer 303. The I FFT transformer 303 has N P outputs 
connected to a selector 305. The selector 305 has a 
plurality of outputs, wherein each succeeding group of Q 
outputs is connected to a corresponding filter and/ ox- 
threshold means 307. Each of the filter means 307 has Q 
outputs connected to a corresponding zero padding means 
309, wherein each of the zero padding means 309 has N c 
outputs connected to a corresponding N c -point FFT 
transformer 311 of a plurality of FFT transformers. Each 
of the FFT transformer 311 has N c outputs connected to a 
corresponding inverse window means 313, each of the inverse 
window means 313 having N c outputs. The plurality of the 
filter means 307, the plurality of the zero padding means 
309, the plurality of the FFT means 311 and the plurality 
of the inverse window means 313 are part of a means 315 for 
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5 determining channel information of the plurality of 

channels based on the selected groups in accordance with a 
further embodiment of the present . invention. 

The N P discrete values of the input signal are first 
10 windowed by the frequency domain window means 301 in order 
to avoid leakage affects which can be caused by the 
following IF FT transformation performed by the IFFT 
transformer 303. The inventive apparatus comprises single 
means 301 for windowing the input signal to obtain a 
15 windowed input signal to be used as the input signal for 

trie transformer. The input signal includes a plurality of 
pilot sequences having different phase shifts with respect 
to each other, wherein each of the pilot sequences 
corresponds to a transmitting point, transmitting the pilot 
20 sequence for channel estimation. 

The IFFT transformer 303 translates the phase shifts into 
delays by performing an IFFT transform applied to the N P 
values of the input signal in order to obtain a transformed 
2 5 signal. It is to be remarked here that the IFFT 

transformer 3 0 3 performs a single IFFT algorithm applied to 
the input signal. Hence, che transformed signal consists 
of successive N P discrete values being related to the 
delays . 

30 

The selector 305 selects a plurality of groups from a" set 
of discrete values of the transformed signal, wherein each 
selected group of the selected groups actually includes 
channel information in the form of a channel impulse 
35 response of each communication channel. 
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5 Due to a possible channel noise , the channel information is 
corrupted by a distortion signal related to the channel 
noise. Kence, the selected .groups are filtered by the 
respective filter means 307 in order to estimate the 
channel information from the corrupted channel information. 

10 The filter means 307 can be operative to perform a 

detection of the channel information from the noisy channel 
information. This can be performed by applying a minimum 
mean squared estimation by filtering the selected groups, 
wherein each filter means 307 comprises a filter having 

15 filter coefficients , wherein each filter coefficient can be 
determined from a multiplication of a first matrix by a 
second matrix, the first matrix being a cross-correlation 
matrix between channel information and the corrupted 
channel information and the second matrix, being an inverse 

20 of an auto-correlation matrix of the corrupted channel 

information as described above. Alternatively,, a simple 
threshold operation can be performed to estimate the 
channel information by means of hard detection. Hence, the 
filter means 307 outputs a plurality of channel impulse 

25 responses as channel information corresponding to each 
communication cnannel . In order to obtain a channel 
transfer function of a particular channel impulse response, 
the means 315 for determining channel information comprises 
means for transforming the channel impulse responses into 

30 frequency domain. 

In its basic form, DFT based interpolation is performed 
simply by adding N c ~Q zeros to the channel impulse 
response estimate, thus extending the length of each 
35 channel impulse response to H c samples by means of zero 
padding by for example attaching the N c -Q zeros at the end 
of each channel impulse response . Hence, after zero padding 
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5 performed by the respective zero padding means 309, a 
length of each channel impulse response is extended to a 
total length of N c . In the next step, the extended channel 
impulse responses in the time domain are transformed into 
frequency domain by applying for example a fast Fourier 

10 transform to each extended channel impulse response by the 
FFT means 311 . In order to obtain a channel transfer 
function corresponding to each channel impulse response, an 
inverse window operation is performed in order to remove 
the windowing effects caused by the frequency domain window 

15 means 301 . 

Fig. 4 shows a conceptual block diagram of an inventive 
apparatus for providing a. plurality of different pilot 
sequences to be used for channel estimation. 

20 

The apparatus shown in Fig. 4 comprises means 401 
calculating an i th value of a piloL sequence to be 
transmitted by a u th transmitting point of N T transmitting 
points , 

25 

In order to calculate distinct, pilot sequences the means 
401 uses the following equation for calculating an i th 
value J£f M) : 

3 0 X iM) = e -;M^~l)^r Jj-'M 

wherein N T denotes a number of transmitting points being 
greater than or equal to 3 and wherein Xf ] denotes an i th 
value of a scrambling sequence, wherein a length of the 
35 scrambling sequence is equal to a length of the pilot 
sequence . 
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5 Hence, the pilot sequence corresponding to a first 
transmitting point, wherein u=l, is identical to the 
scrambling sequence. A n th pilot sequence u=n to be 
transmitted by an n tn transmitting point comprises a number 
of discrete values related to the phase shifts of the n th 

10 pilot sequence,., wherein each value of the n th pilot 
sequence is multiplied by the scrambling sequence. The 
phase shifts of each value of the pilot sequence depend on 
a number of the value being calculated, wherein the number 
of the value is associated with a position of the value 

15 within the pilot sequence. Furthermore, the phase shift of 
the n th pilot sequence with respect to the other pilot 
sequence depends on a. number of transmitting points divided, 
by a total number of transmitting points N T . Hence, 
succeeding values of each pilot sequence have a phase shift 

20 with respect to each other, wherein the phase shift is 
negative in a mathematical sense, due to a negative sign 
appearing in the exponent of the exponential function. 

Hence, the pilot sequences provided by the inventive 
2 5 apparatus are suited for low complexity DFT based channel 
estimation. By using the inventive phase shifted pilot 
sequences, the properties of the Fourier transform can 
efficiently be exploited to separate the N T superimposed 
signals. It is well known that a phase shifted signal in 
30 the frequency domain is translated to a delay in the time 
domain by the Fourier transform and vice versa. So,' by 
using phase shifted pilot sequences in the frequency 
ccmain, the resulting time domain signal will be delayed. 
Furthermore, by assigning a certain phase shift as 
35 described above to each transmit antenna' s pilot, the 
corresponding delay in the time domain will be dependent on 
the transmit antenna. Provided that a length of each 
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5 training sequence is greater than or equal to K T Q f whore Q 
denotes a length of a channel having a maximum length of 
the plurality of the communication channels to be 
estimated, the N T superimposed signals can be separated 
completely simply by performing an inverse discrete Fourier 
10 transform. 

It has been discussed above that phase shifted pilot 
sequences provided by the inventive apparatus shown in Fig, 
4 are suited to time domain MIMO OFDM channel estimation 
15 since these sequences are the perfect match to the IDFT 
operation. Nevertheless, other orthogonal sequences may 
also be applied to a low complexity channel estimation with 
some additional processing in time domain. 

20 Fig. 5 shows a conceptual block diagram of a further 
inventive apparatus for providing a number of different 
pilot sequences to be used for channel estimation, wherein 
the pilot sequences are orthogonal sequences based on 
columns of a Hadamard matrix, 

25 

The apparatus shown in Fig. 5 comprises a means 501 for 
generating a Hadamard matrix having a number of columns f 
wherein the number of columns corresponds to rhe number of 
different pilot sequences (pilot symbols) . The means 501 

30 has a. number of outputs,, the number of outputs being equal 
to the number of columns of the Hadamard matrix. " The 
outputs of the means 501 are connected to a means 503 for 
calculating different pilot sequences. The means 503 for 
calculating different pilot, sequences has a number of 

35 outputs for providing the different pilot sequences. 
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5 The means 501 for generating a Hadamard matrix calculates 
columns of the Hadamard matrix being orthogonal to each 
other . The number of columns of the Hadamard matrix is 
provided to the means 503 for calculating different pilot 
sequences, wherein the means 503 calculates the different 

10 pilot sequences on the basis of the columns of the Hadamard 
matrix. Each pilot sequence is derived from, another column 
of "he Hadamard matrix. In order to generate a pilot 
sequence, the corresponding column of basis functions given 
by the columns of the Hadamard matrix can for example be 

15 periodically extended to a pilot sequence length. 

However, the inventive pilot sequences based on Hadamard 
codes cannot be applied to the simplified receiver 
structure discussed above in a straight forward way for a 

20 number of transmitting points being greater than 2. With 
some additional processing a channel estimation scheme 
based on Hadamard codes can be simplified. For example, if 
four transmitting points are applied to transmit pilot 
sequences, then a channel information determined by a DFT 

25 based channel estimation scheme has an additional phase 
shift factor which is a function of a phase shift of 180 
degrees divided by four. Hence, the additional processing 
necessary to determine channel information is small since 
only phase shift operations , additions and subtractions are 

30 required in order to solve a set of linear equations. 

If, for example, eight transmitting points are applied to 
transmit pilot sequences, the estimated channel information 
has an additional phase shift factor which is a function of 
35 a phase shift by 180 degrees divided by the total number of 
transmitting points which is eight in this case. Hence, 
only four phase shifting and four add and subtract 
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5 operations are necessary. Note that application to Walsh 
sequences is straight forward since Walsh sequences are 
obtained by rearranging the columns of the Hadamard matrix, 
i.e. assigning different pilot sequences to the transmit 
antennas . 

10 

Pilot aided channel estimation (PACE) was first introduced 
for single carrier systems and required a flat fading 
channel, as described by J . K . Cavers, "An Analysis of Pilot 
Symbol Assisted Modulation for Rayleigh Fading Channels," 

15 IEEE Transactions on Vehicular Technology, vol. VT-40, pp. 
686-693, Nov. 1991. It is based on periodically inserting 
known symbols, termed pilot symbols in the data sequence. 
If the spacing of the pilots is sufficiently close, to 
satisfy the sampling theorem, the channel estimation and 

20 interpolation for the entire data sequence is possible. 
When extending the idea of PACE to multi-carrier systems, 
it must be taken into account that for OFDM fading 
fluctuations in two dimensions occur, i.e., in time and 
frequency. In order to satisfy the two-dimensional sampling 

25 theorem,, the pilot symbols are scattered throughout the 
time/frequency grid, which yields to a two-dimensional 
pilot grid. 

Possible realizations of pilot grid structures for OFDM 
30 channel estimation are illustrated in Fig. 6. One 
possibility to structure a pilot, grid is to transmit ' one 
OFDM symbol containing only pilots followed by D t -1 data 
symbols. This sequence is shown in the left diagram in Fig. 
6. This scheme is applicable for channels with little time 
35 variations, as observed in an indoor environment. In this 
case, no interpolation (in frequency direction) is 
necessary. Another possibility is to transmit the pilots on 
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5 reserved carriers as shown in the middle diagram in Fig. 6. 
This scheme can support mobility but requires interpolation 
in frequency direction. A .more efficient solution is to 
employ a scattered pilot grid, as shown in the right 
diagram in Fig. 6. The structure for OFDM signalling allows 

10 a channel estimator to use both time and frequency 
correlations. With such a scattered pilot grid, the 
property of OFDM, thac adjacent sub-carriers as well as 
adjacent symbols are correlated, can be exploited for 
channel estimation. Thus, the overhead due to pilots can be 

15 reduced but interpolation in time and frequency is 
required. 

For pilot symbol aided channel estimation (PACE) based on a 
scattered pilot grid, known pilot symbols are multiplexed 

20 into the data stream. Interpolation is used to obtain the 
channel estimate for the information carrying symbols. PACE 
for single carrier systems was introduced by Cavers in the 
above mentioned prior art document. P. Hoeher, S. Kaiser, 
and P. Robertson, "Two-Dimensional Pilot-Symbol-- Aided 

25 Channel Estimation by Wiener Filtering", in Prcc. IEEE 
Intern. Conf. On Acoustics, Speech, and Signal Processing. 
(ICASSP'97), Munich, Germany pp. 1845-1848, 1997, and R. 
Nllsson, 0. Edfors, M. Sandell, and P, Borjesson, "An 
Analysis of Two-Dimensicnal Pilot-Symbol Assisted 

30 Modulation for OFDM", in ?roc. IEEE Intern. Conf. On 
Personal Wireless Communications (ICPWC 97) , "wumbai 
(Bombay) , India, pp. 71-74, 1997 . They propose two- 
dimensional 2D filtering algorithms for PACE. Such a 2D 
estimator structure is generally too complex for practical 

35 Implementation. To reduce the complexity, separating the 
use of time and frequency correlation has been proposed. 
This combined scheme, termed double 1-dimensiona.l PACE uses 
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separate Wiener filters f one in frequency direction and one 
in time direction. 

A general formal description of a regular grid in a 2D 
plane can be found in Y. Li, "Pi 1 o t - S ymbc 1 - Ai d e d Channel 
Estimation for OFDM in Wireless Systems", IEEE Transactions 
on Vehicular Technology,, vol. 49, pp. 12C7-I215, July 2000. 

The total number of pilots transmitted by one frame is 
depicted as N' p ~N c f D f and N p = L/D ff representing che 
number of pilots in frequency and time direction, 
respectively. The following notation will be used. Given a 
variable describing a 2D structure X, the subsets which 
describe the dimension corresponding to frequency and time 
directions are denoted by X' and X' ' , respectively. As a 
general contention, variables describing pilot symbols will, 
be marked with a in. the following 

For MIMO systems every transmit antenna signal uses its own 
pilot, to enable the receiver to separate the superimposed 
signals. To describe pilot symbol assisted channel 
estimation, it is useful to define a subset of the received 
signal sequence containing only pilots, jx£'l!j = {xjf} with 
l = lD t and i = i D f . So, the pilot sequences transmitted at D f 
times lower rate i = \j!D f \ in frequency direction, and a D t 
times lower rate / =[//£),] in time direction, respectively. 
It is assumed that the pilots \x$\ are chosen from a phase 
shift keying (PSK) constellation |jf^j=l. Before 
transmission, the pilots \x-fj\ are multiplied by an outer 
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5 pilot sequence |? 0j „j which is identical for all transmit 
antennas tc yield the transmitted pilot sequence . 




10 The outer pilot sequence is chosen to have a low peak 

to average power ratio (PAPR) in the time domain and/or to 
have good correlation properties for synchronisation etc. 
In any case, the outer pilot sequence is introduced solely 
for the sake of completeness. 

15 

At the receiver, the cyclic prefix is removed and an FFT is 
performed to yield the received signal after OFDM 
demodulation. Assuming perfect synchronisation, the 
received signal Y u is obtained. For channel estimation, 
20 the received signals at the pilot positions are 
demultiplexed from the data stream, and after removing the 
outer pilot sequence , by dividing through X hy the received 
pilot is obtained. 

2 5 Y TJ = F WeG -■- Xl£ G H}£ s + N !j€G 

where G is the subset of the OFDM frame containing the 
pilots. 

30 Considering the pilot sequence of OFDM symbol 1 = Td, from 
transmit antenna u which can be expressed by a column size 
factor of size Np 
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where the transmitted pilot sequence, the channel transfer 
function (CTS) and additive r.cise terra, are given by 

10 



15 

r.he N' p x A r ^ D FT -matrix F transforms the CIR into the 
frequency domain, defined by 

20 ^e i2m ' iN >; Q<i<N' P 1, Q<n^N' P -l 

In case that Q<N' F , the last N' P ~Q of the DFT output need 
to be removed, which can be formally performed by the 
matrix I y ^ e - [^^Qfi.^-Q^gt of dimension N' P x Q , with entries 
2 5 equal to 1 at the main diagonal and 0 elsewhere. Fcr 
Q~y p , the matrix i.^^ becomes the identity matrix. In 
practice, the DFT transformation can be efficiently 
generated using a N' p -point FFT. In case that Q<N' pf the 
last N'p -Q outputs are skipped. 
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5 

Hence, the flowing equation follows 
Y{ - X'jHL + Nl = X' 7 F Nt K + Nj 
10 where 



15 



Referring now tc a design of phase shifted pilot sequences, 
20 the set of N T phase shifted pilot sequences of length N' p 
is defined by 

jW^^MHH)/^. ^ = {l,.- v V r } 5 T = il,-,N>} 

25 wherein the pilot sequences defined by the above equation 
are in accordance with the preferred embodiment according 
to the present invention . 

In stead of pre-multiplying the received signal prior to 
30 the ID FT operation, the received pilot signal of OFDM 
symbol l = lD ft Y~ is directly transformed into the time 
domain : 
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5 

r Li p"?; = ±-Yf f! x^F! M . +±~F H N' eC"'' 1 

Hence, only one IDFT is required to generate £t 

10 As to the properties of the DFT, phase shifted sequences 
are transformed by the DFT to a time shift. By assigning 
each transmit antenna its own phase shift, the time dona in 
signal g 7 is separated perfectly if Q<N' P /N r . 

15 In order to analyse the following equation, it is useful to 
examine the components of the vector £ r appearing in the 
above equation. The n~ h entry of c,' f is in the form 



25 
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Where n 7 g is the IDsl of N TT . Hence, calculation of the 

DFT of the phase shifted pilot sequences in accordance with 
the preferred embodiment of 'the present invention leads to- 

It is noted that this is identical to the DFT of X?\*Xj% . 
since Xfl * is the all one sequence. Hence, the desired 
result for £ T ^ is obtained by 

with n = q \ (u~r)N' F !N T and Q<N' P !N T 

If we assume that Q = N' P /N T , column n = q + (ju-1)Q of the CIR 
is given by Otherwise, if Q<N' P iN Tr the size of 

the channel impulse response can be matched by zero 
padding. Hence , the following equation follows readily 

4>^ + «^/4 T n;>Av2 

where n[ =\i N' P F" N' { denotes the I DFT of the AWGN realisation 
N'j , which is itself a white Gaussian noise process with 
the same variance. This means that the N T superimposed 
signals are completely separated, provided that N' P >N T Q . 
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The MMSE estimator is in general dependent on the choice of 
the pilot symbols. However, choosing appropriate pilot 
sequences, the estimator _ becomes independent of the 
transmitted pilots. It is desirable to choose a set of 
pilot sequences , which minimises the mean squared error 
(i.e. the performance of the estimator} and. the 
computational complexity of the estimator. Hence, if the 
pilot sequences are appropriately chosen, the least squares 
estimator as well as the minimum mean square error 
estimator can be grossly simplified. 

In order to generate the MMSE estimator, a .knowledge of the 
correlation matrices R'., and R' hc are required. Fox the 
derivation of the Wiener filter it is assumed -hat 
Q-N' P l N T . Otherwise, the size of the channel impulse 
response can be matched by zero padding. The MMSE estimates 
for OFDK symbols is given by (S.M, Kay, "Fundamentals of 
Statistical Signal Processing: Estimation Theory", 
Englewood Cliffs, NJ; Prentice Ea.ll International, 1393) 



w' = R' h( R£ 



where the correlation matrices R' (i and R' h .- are defined by 
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and 

10 
15 

In order to perform MIMO OFDM channel estimation in cime. 
direction, by appropriately scaling of the parameters, the 
same inventive algorithms for channel estimation can be 
applied to D FT -based channel estimation in time -direction. 

20 

In the following, the system model for equalisation of N T 
superimposed pilot signals in time direction, i.e. the 
received pilot sequence of tone 7 of N'p OFDM symbols is 
introduced. Considering channel estimation in time 
25 direction, the pilots [y, ^,7^ r^^] sG are used in. order 
to estimate H\f . In vector notation, the received pilot 
sequence of tone i is defined by 



30 
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where the transmitted pilots, the general transfer 
function, the general impulse response and the additive 
noise components „ are given by 



There are two scenarios which will be considered: first, 
burst like transmission where a frame of an OFDM symbol is 
transmitted; second, a broadcasting scenario with a 
continuously transmitted data stream. In burst transmission 
mode, L OFDM symbols are received,- buffered,- and after 
reception is complete f the whole frame is processed. This 
is the same as for the channel estimation in frequency 
direction, where all pilots of one OFDM symbol are being 
received together as well. The difference is that the 
algorithms which were applied dependent on the i variable 
are now used in the 1 variable. That is, the received pilot 
is transformed into the Doppler domain via an IFFT . The 
least squares and minimum, mean squared error estimator work 
exactly in the same way as the previously discussed case of 
channel estimation in frequency direction. 
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5 Referring to the broadcasting scenario, for continuous 
transmission there are some differences with respect to the 
channel estimation in frequency direction. Some buffering 
is usually possible since receiving the whole sequence 
before processing is soir.e times not feasible. L = D t N'p is 

10 defined as the number of OFDM symbols, which are in the 
input buffer. For smoothing type filtering as described by 
S.M. Ray, "Fundamentals of Statistical Signal Processing; 
Estimation Theory", Englewood Cliffs, NJ: Prentice Hall 
International, 1993, we use past as well as future pilots 

15 to estimate H\f . This means 1 < / < D t N" . Clearly, for 
smoothing an estimate cannot be obtained until all pilots 
have been received, which requires buffering of A/ =- D s Ny - 1 
OFDM symbols. The alternative is to use prediction type 
filtering where l>D t N'p. In this case only past pilots are 

20 used for channel estimation in time direction. Obviously, 
prediction type filtering does not require any buffering, 
however, the performance with respect to smoothing 
degrades , 

2 5 The major difference with respect to channel estimation in 
frequency direction is that only D t symbols per sub-carrier 
are output at a time, i.e., an input of Np pilots is used 
to produce an output of D t symbols per sub-carrier . 
Furthermore, the number of pilots in time direction is 

30 usually significantly smaller than the number of pilots in 
frequency direction, but the pilot spacing frequency in 
time direction is similar. Since Np«N' F the leakage of 
the Doppler domain transform will be more significant. On 
the other hand, when only the symbols in the centre of the 

35 processed pilots are used, the performance may be improved, 
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since the estimation error is typically much lower at the 
centre compared to the area near the edges . 

Up to now much of discussion is dealing with ID channel 

estimation in frequency direction over one OFDM symbol in 
the i variable. .The proposed techniques are, however, also 
applicable to channel estimation in time direction in a 
straight forward way (i.e. estimating in the 1 variable). 
Moreover, an extension to a multi carrier system with 
correlations in two dimensions is possible by using two 
cascaded one dimensional Wiener filters, This idea has beer. 
suggested, for SISC channel estimation by e.g. P. Hoeher, S . 
Kaiser, and P. Robertson, "Two-Dimensional Pilot "-Symbol- 
Aided Channel Estimation by Wiener Filtering", in Proc 
IEEE Intern. Conf . On Acoustics, Speech, and Signal 
Processing. (ICASSP' 97) , Munich, Germany pp. 1845-1848, 
1997, and can be extended to the MIMC case. 

The inventive apparatus for estimating a plurality of 
channels is operative to translate phase shifts into delays 
in order to estimate channel information, i.e. a channel 
impulse response of one of a plurality of communication 
channels . Hence, phase shifted sequences are perfectly 
suitable for channel estimation. Moreover, also other 
orthogonal sequences can be applied to channel estimation, 
as long as the pilot symbols have a particular phase shift 
characteristic with respect to other pilot symbols . That 
means, that the phase shift between two pilot sequences may 
vary following a. predetermined function so that for example 
a first half of a first pilot sequence has a phase shift 
with respect to a first half of a second pilot sequence and 
a second half of the pilot sequence has a further phase 
shift with respect to a second half of the second pilot 
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sequence . Therefore, -he phase shift characteristics of the 
pilot sequences can be exploited at the transmitting points 
to estimate the channel -by- applying for example the 
inventive low complexity channel estimation scheme. 

As an example o.f orthogonal sequences, the properties of 
Haaamard sequences are studied. For a system with N T 
transmit antennas a set Hadamard basis function of order 
N T , denoted by H is used. In order to generate X'-^ , the \x tk 
column of these basis functions are periodically extended 
-co pilot sequence length N' pr as discussed above. For 
example, for N~-2 the Hadamard sequences are identical to 
the phase shifted sequences, i.e. Jif~ (l) is the all one 
sequence, and the second sequence is oscillating between 1 
and -1, i.e. X f } =[l-l,l,--lf . If N T »4 r the Hadamard 
sequences can be cast into the matrix 



1111 

1-1 1-1 
1 1 -1 -1 
1 -1 -1 1 



It is easy to clarify that Hi is symmetric and orthogonal 
i.e. II A = Hi and H 4 -H\ =4-2 . In case that N T =4 f the received 
pilot sequence transformed into the time domain, is defined 
by £7 = \l NpF"Y~ . It will be shown, that the least squares 
estimators for transmit antennas 1, 4 is given by 
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fey - e^Xf ) + e- J " ! X? = Tf ] + HI 



10 where Cj' } accounts the part, of £ T , assuming that Is 

divided into N T equal parts, i.e. £ ? = \l f t9 '} . 

Furthermore , denotes the resulting noise term of u th 

part. It is seen that the signal of antenna 2 is received 
in partition 3. This may be confusing and can be fixed by 
15 exchanging columns 2 and 3 of H NT . 

To apply the above simplification of the receiver structure 
with pilot sequences based on Hadamard codes, additional 
processing may be necessary. For the case N T =4 and pilot 

2 0 sequence is based on Hadamard codes, the receiver structure 
of the FFT Is derived. The starting point Is the received 
pilot sequence transformed into time domain given by 
£ J =liN , p F H Y{ , The DFT of the pilot sequence Xf^ for 
transmit antenna 1, 4 can 

25 be expressed as 
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10 Hence, it can be observed that the signal from antenna I is 
identical to the case of phase shifted sequences are used, 
and the signal from antenna 3 is the same as the signal 
from antenna 2 in the above equation, which corresponds to 
the least squares (LS) estimate. This is due to the 

15 inventive arrangements of the codes and may be fixed by- 
exchanging for example columns 2 and 3. The DFT of pilots 3 
and 4 in the -above equation consists of two impulses which 
interfere with each other. Furthermore, the DFT of pilot 3, 
F ff'Mfj\^f is the conjugate complex of F K (xfl) . Hence, 
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5 where again it is assumed that Q < N? I N r . Removing the 
inserted phase shift of the pilots and adding the two-parts 
of the signal, the equation to- be- shown follows readily. 

In practice, two times oversampling provides a good trade 
10 of between minimising the system overhead due to pilots and 
optimising the performance. It is assumed that the length 
of the guard interval is not smaller than Q, i.e. the guard 
interval is not . shorter than the maximum delay of the 
channel . 

15 

Moreover, depending on certain implementation requirements 
of the inventive method for estimating a plurality of 
channels or for providing a plurality of different pilot 
sequences can be implemented in hardware cr in software. 

20 The implementation can be performed using a digital storage 
medium, in particular a disc or a CD having electronically 
readable control signals stored thereon, which can 
cooperate with a programmable computer system such that the 
inventive methods are performed. Generally, the present 

25 invention is, therefore, a computer-program product with a 
program code stored on a machine- readable carrier, the 
program code being for performing the inventive methods, 
when the computer program product runs on a computer . In 
other words, the inventive method is, therefore, a computer 

30 program having a program code for performing the inventive 
methods, when the computer program runs or. a computer. 
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5 CLAIMS 



1. Apparatus for estimating ■ a plurality of channels 
from a plurality of distinct transmitting points to 
a receiving point, wherein an input signal includes 

10 a superposition of a first signal based on a first 

pilar sequence transmittabie from a first 
transmitting point and a second signal based on a 
second pilot sequence transmittabie from a second, 
transmitting point, the first pilot sequence and 

15 the second pilot sequence being orthogonal to each 

other within, a predetermined orthogonality range 
and having a phase shift to each other, comprising: 

a transformer (103; 303) for transforming the input 
signal into a transformed signal having a set of 
discrete values, the transformer (103; 303) being- 
operative to translate the phase shift between the 
first and the second pilot sequence into a delay, 
so that a first group of the set of the discrete 
values includes channel information of a first 
channel from the first transmitting point to the 
receiving point and that a second group of the set 
of discrete values includes channel information of 
a second channel from the second transmitting point 
to the receiving point ; 

a selector (105; 305} for selecting the first, group 
from the set of discrete values and for selecting 
the second group from the set of discrete values to 
35 obtain selected groups; and 

means (107; 315) for determining channel 



25 
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5 information for the plurality of channels based on. 

the selected groups . 

2. Apparatus in accordance with claim 1, wherein the 
input signal includes,, in a superimposed manner, a 

10 third signal based on a third pilot sequence 

t ran smitt able from a third transmitting point, the 
third pilot sequence being orthogonal to the first 
and the second pilot sequences within the 
orthogonality range, and having a further phase 

15 shift to the first pilot sequence, the transformer 

(103; 303) being further operative to translate the 
further phase shift into a further delay, so that a 
third group of the set of discrete values includes 
channel information of a third channel from the 

20 third transmitting point to the receiving point ; 

and 

wherein the selector (105; 305) is operative to 
further select the third group from the set of 
25 discrete values to obtain a selected third group. 

3. Apparatus in accordance with claim. 1 or 2, wherein, 
each group consists of an equal number of 
successive discrete values of the set of discrete 

30 values . 

4. Apparatus in accordance with any one of the 
preceding claims , wherein the transformer is 
operative to perform a Fourier transform or a fast 

35 Fourier transform or a discrete Fourier transform 

or an inverse Fourier transform or an inverse fast 
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5 Fourier transform or an inverse discrete Fourier 

transform on the input signal . 

5. Apparatus in accordance with any one of the 
preceding claims, wherein the set of discrete 

10 values includes a total number of discrete values, 

the number being obtained by the transformer (103; 
303) . 

6. Apparatus in accordance with claim 5, in which the 
15 plurality of transmitting points has a total number 

of transmitting points, wherein the selector (105; 
305} is operative to select, as a group, a number 
of successive discrete values, the number being 
equal to the total number of discrete values 
20 divided by the total number of transmitting points. 

7. Apparatus in accordance with claim 6, in which the 
selector (105; 3C 5 ) is operative to select the 
groups such that a first discrete value of the set 

25 of discrete values is a first discrete value in the 

first group and a first discrete value of uhe 
second group of discrete values is a discrete value 
in the set of discrete values succeeding a discrete 
value having a highest order of the values selected 

30 to obtain the first group. 

8. Apparatus in accordance with any one of the 
preceding claims, in which the input signal is a 
set of successive input signal values, and in. which 

35 the transformer (103; 303) is operative to 

transform the set of successive input signal values 
into the set of discrete values, wherein a total 
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number of discrete values of the set of discrete 
values is equal to a total number of input signal 
values in the set of. input signal values. 

9. Apparatus in accordance with any one of the 
preceding claims, in which the transformer (103; 
303) is operative to perform, a transform algorithm, 
and in which the set of discrete values is 
generated by a single execution of the transform 
algorithm by the transformer. 

10 . Apparatus in accordance with claim 9 , in which the 
transform algorithm being performed by the 
transformer (103; 303) is a Courier transform or a 
fast Fourier transform or a discrete Fourier 
transform or an inverse Fourier transform or an 
inverse fast Fourier transform, or an inverse 
discrete Fourier transform algorithm applied to the 
total number of input signal values . 

11. Apparatus in accordance with any one of the 
preceding claims, wherein the pilot sequences 
include a complex amplitude sequence and phase 
information, wherein the complex amplitude sequence 
is the same for the pilot sequences, comprising; 

a multiplier (101) being operative to multiply a 
received signal by a complex-conjugated version of 
the amplitude sequence to obtain the input signal. 

12. Apparatus in accordance with claim 11, wherein the 
received signal is a set of successive discrete 
received signal values and wherein the complex 
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5 conjugated version of the amplitude sequence is a 

set of successive discrete complex conjugated 
amplitude signal values., and 

wherein the multiplier (101.) is operative to 
10 perform a single complex multiplication for a value 

cf the set of successive discrete signal values and 
a corresponding value of the set of successive 
discrete complex conjugated amplitude signal 
values, so that after a number of single 
15 multiplications the set of successive discrete 

input values is obtained. 

13. Apparatus in accordance with any one of the 
preceding claims, in which the means (107; 315) for 

2 0 determining channel information is operative to 

receive the selected groups and to output the 
selected groups, each group of the selected groups 
including channel information of a particular 
channel . 

25 

14. Apparatus in accordance with claim 13, wherein the 
means (107; 315) for determining channel 
information is operative to determine the channel 
information by output ting the selected groups,, each 

3C group including channel impulse responses of a 

particular channel as channel information. 

15. Apparatus in accordance with claim 13, wherein the 
means (107; 315} for determining channel 

35 information includes a time/frequency conversion 

means for obtaining a channel transfer function as 
the channel information. 
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16. Apparatus in accordance with any one of the 
preceding claims, wherein- a group of the selected 
groups and a further group of the selected groups 
include channel information of a channel and 

.0 channel information of a further channel in a super 

imposed manner, wherein the means (107; 315) for 
determining channel information is operative to 
determine the channel information of the channel 
and the channel information of the further channel 

.5 based on the group and on the further group. 

17. Apparatus in accordance with claim 16, wherein a 
group of the selected groups includes a 
superposition of the channel information of a 

:0 channel being multiplied by a phase shift factor 

and of a further channel information of a further 
channel being multiplied by a further phase shift 
factor, and 

!5 wherein a further group of the selected groups 

includes a superposition of the channel information 
of the channel being multiplied by the complex 
conjugated phase shift factor and of the further 
channel information of the further channel being 

JQ multiplied by the complex-conjugated further phase 

shift factor, the phase shift factor arid ' the 
further phase shift factor being caused by the 
pilot sequences, wherein 



35 



the means (107 ; 315) for determining channel 
information is operative to determine the channel 
information of z.he channel and the channel 
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information of the further channel using the group 
and the further group. 

18. Apparatus in accordance with claim 17, wherein the 
means (1C7; 315) for determining channel 
information is operative to determine the channel 
information and the further channel information 
using phase shift operations and addition and /or 
subtraction operations applied to the group and to 
the further group. 

19. Apparatus in accordance with claim 17 or 18, 
wherein the phase shift factor is a function of a 

phase shift of 180 degrees divided by the total 
number of transmitting points. 

20. Apparatus in accordance with any one of the 
preceding claims, wherein a. channel information is 
corrupted by a distortion signal, the means (107; 
315) for determining further comprising an 
estimator being operative to estimate the channel 
information from the corrupted channel information. 

21. Apparatus in accordance with claim 20, wherein the 
estimator is operative to perform a minimum mean 
squared error estimation or a maximum a posteriori 
estimation or a maximum likelihood estimation for 
providing an estimated channel information. 

22. Apparatus in accordance with any one of the 
preceding claims, in which the pilot sequences are 
used to modulate a discrete carrier in a multi 



WO 2IMI4,'HX<i7(>6 



63 



PCT/EP2003/00J211 



5 carrier modulation scheme at a predetermined time 

instant, the apparatus further comprising: 

a controller for detecting the input signal in a 
received signal and for enabling the transformer 
10 1103; 303) , the selector (105; 305) or the means 

(107; 315} for determining only when the input 
signal is detected in the received signal. 

23. Apparatus in accordance with any one of the 
15 preceding claims, in which the pilot sequences are 

used to modulate the same pilot carrier in a multi 
carrier modulation scheme at discinct transmitting 
time instants, the apparatus further comprising: 

20 a controller for detecting a pilot carrier in a 

received signal including successive multi carrier 
modulated symbols and for collecting values of the 
pilot carrier at time instants corresponding tc 
transmitting instants to obtain the input signal,. 

25 wherein the channel information obtained by 

processing the input signal, using the transformer 
(103; 303} , the selector (105; 305) and the means 
(107; 315) for determining includes information on 
a multi carrier channel from transmitting points to 

30 receiving point defined by the pilot carrier over- 

time . 



24. Apparatus in accordance with any one of the 
preceding claims, wherein the apparatus further 
comprises a single means (301) for windowing the 
input signal to obtain a windowed input signal to 
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5 be used as the input signal for the transformer 

(103; 305) . 

25. Apparatus for providing a plurality of different 
pilot sequences to be used for channel estimation, 
10 comprising: 

means (401) for calculating an 1 th value Xi (p) of a 
pilot sequence to be transmitted by a u th 
transmitting point of M T transmitting points using 
15 the following equation: 



20 wherein M T denotes a number of transmitting points 

and wherein N T is greater than or equal to 3; and 

wherein denotes an i th , value of a scrambling 

sequence, the length of the scrambling sequence 
25 being equal to the length of the pilot sequence. 

26. Apparatus for providing a number of different pilot 
sequences to be used for channel estimation, 
comprising; 

30 

means (501) for generating a Kadamard matrix having 
a number of columns, the number of columns 
corresponding to the number of different pilot 
sequences; 



means (503) for calculating the different pilot 
sequences based on the columns of the Hadamard 
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5 matrix, wherein each pilot sequence is derived from 

another column of the Hadamard matrix. 

27. Method for estimating a plurality of channels from 
a plurality of distinct transmitting points to a 

10 receiving point, wherein an input signal includes a 

superposition of a first signal based on a first 
pilot sequence transm.itta.ble from a first 
transmitting point and a second signal based on a 
second pilot sequence transmittabie from a second 

15 transmitting point, the first pilot sequence and 

the second pilot sequence being orthogonal to each 
other within a predetermined orthogonality range 
and having a. phase shift to each other, comprising 
the following steps: 

20 

transforming the input signal into a transformed 
signal having a set of discrete values by 
translating the phase shift between the first and 
the second pilot sequence into a delay,, so that a 

2 5 first group of the set of discrete values includes 

channel information of a first channel from the 
first transmitting point to the receiving point and 
that a second group of the set of discrete values 
includes channel information of a second channel 

30 from the second transmitting point to the receiving 

point; 

selecting 7. ho first group from the set of discrete 
values and selecting the second group from the set 
35 of discrete values to obtain selected groups; 
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determining channel information for the plurality 
of channels on the basis of the selected groups. 

Method for providing a plurality of different pilot 
sequences to be used for channel estimation, 
comprising the following steps: 

Calculating an i th value value Xi t;i) of a pilot 
sequence to be transmitted by a u th transmitting 
point of I\V transmitting points using rhe following 
equation : 



wherein N T denotes a number of transmitting points 
and wherein N t is greater than or equal than 3; and 

wherein xf ] denotes an i th value of a scrambling 
sequence, the length of which being equal to the 
length of the pilot sequences . 

29. Method for providing a number of different pilor. 
sequences to be used for channel estimation, 
comprising the following steps: 

generating a Hadamard matrix having a number of 
columns, the number of columns corresponding to the 
number of different pilot sequences; 

calculating the different pilot sequences based on 
the columns of the Hadamard matrix, wherein each 
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pilot sequence is derived from another column of 
the Hadamard matrix. 

30. Computer program having a program code for 
performing the method for estimating a plurality of 
channels., in accordance with claim 27, the method 
for providing a plurality of pilot sequences in 
accordance with claim 28 r the method for providing 
a plurality of pilot sequences in accordance with 
claim 29, when the program runs on a computer. 
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